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Dynamical calculations are performed for all isomultiplets of the flavour antidecuplet 



to which the newly discovered pentaquark G + belongs. The framework is a constituent 
quark model where the short-range interaction has a flavour-spin structure. In this model 
the lowest pentaquarks have positive parity. Each antidecuplet member is described by a 
variational solution with the Pauli principle properly taken into account. By fitting the 
mass of B + of minimal content uudds, the mass of S , of minimal content ddssu, is 
predicted at approximately 1960 MeV. The influence of the octet-antidecuplet mixing on 
the masses of the Y = 1 and pentaquarks is considered within the same model and the 
role of the kinetic energy plus the hyperfine interaction in this mixing is pointed out. 



1. Introduction 
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The existence of exotic baryons containing four quarks and an antiquark in their lowest 
Fock component has now a solid experimental support. The observation of a narrow peak 
at 1.54 ± 0.01 GeV/c 2 , called G + , as an S = 1 baryon resonance in the photo-production 
from neutron — > K + K~n [[Q, has been confirmed by several groups in various photo- 
nuclear reactions [ E] • This has been followed by the observation in pp collisions [ El of 
other narrow resonances H and S° at about 1862 MeV, from which H is interpreted 
as another pure exotic member of an SU(3) flavour antidecuplet. The work of Diakonov, 
Petrov and Polyakov [0] has played a particularly important role in these discoveries. In 
the context of a chiral soliton model they predicted a narrow pentaquark, with a width 
of less than 15 MeV, located at the about experimentally observed mass of 0. 

At the end of the '70's, following the observation of several signals, light pentaquarks 
were studied theoretically [ El E] , but these signals were not confirmed. Charmed pen- 
taquarks with strangeness, uudsc and uddsc were also predicted [E|E], but experimental 
searches carried at Fermilab have remained inconclusive [ E] • These pentaquarks were 
introduced in the context of the one-gluon exchange model (colour-spin interaction) and 
the heavy ones carried negative parity. On the other hand positive parity pentaquarks 
containing heavy flavours were proposed in the context of a pseudoscalar exchange model 
(flavour-spin interaction) [E] about ten years later [H21- in this model, the lowest ones, 
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uuddc and uuddb, do not carry strangeness. Recently the HI Collaboration at DESY [ 
.10] reported a narrow resonance of mass 3099 MeV, interpreted as a uuddc pentaquark. 

The spins and parities of G + and H are not yet known experimentally. In this new 
wave of pentaquark research, most theoretical papers take the spin equal to 1/2. The 
parity is more controversial. In chiral soliton or Skyrme models the parity is positive 
[Hj. In constituent quark models it is usually positive. In the present approach, the 
parity of the pentaquark is given by P = (— ) + 1 , where £ is the angular momentum 
associated with the relative coordinates of the g 4 subsystem. We analyze the case where 
the subsystem of four light quarks is in a state of orbital symmetry [31] and carries an 
angular momentum i = 1. Although the kinetic energy of such a state is higher than that 
of the totally symmetric [4] G state, the [31] G symmetry is the most favourable both for 
the flavour-spin interaction [E| and the colour-spin interaction [E]- in the first case the 
statement is confirmed by the comparison between the realistic calculations for positive 
parity [E] and negative parity [E], based on the same quark model [E]- I n Ref- [E] 
the antiquark was heavy, c or b, and accordingly the interaction between light quarks and 
the heavy antiquark was neglected, consistent with the heavy quark limit. In Ref. [ ITE] 
an attractive spin-spin interaction between s and the light quarks was incorporated and 
shown that a stable or narrow positive parity uudds pentaquark can be accommodated 
within such a model. This interaction has a form that corresponds to rj meson exchange 
[ El and its role is to lower the energy of the whole system. 

The purpose of this letter is to perform dynamical calculations of all the members of 
the antidecuplet to which + and S are supposed to belong. To our knowledge this is 
the first attempt in this direction. The present study is a natural extension of Ref. [ E] 
where the heavy antiquark c or b is now replaced by a light quark u, d or s. To describe 
the short range interaction we rely on the same model [El as that used in [E]- That 
means that the quark-quark interaction has a flavour-spin structure [ E] and that the 
parameters are fitted to the light non-strange and strange baryon spectra. Moreover we 
assume that the quark- antiquark interaction is proportional to a spin-dependent operator, 
but it is flavour independent, as in Ref. [El- Its role is to introduce the same flavour 
independent shift for each member of the pentaquark antidecuplet of equal spin. We shall 
fix this shift by adjusting the mass of + to the experimental value. There is no other 
free parameter in the Hamiltonian model used in this study. For the pure exotic H , 
we predict a mass of 1960 MeV. For the antidecuplet members with Y = 1 and we 
investigate the role of the octet-antidecuplet mixing. To some extent this study will be a 
comparative one. 

We search for a variational solution of a five-body Hamiltonian, containing a kinetic en- 
ergy term, a confinement term and a short range (hyperfine) interaction having a flavour- 
spin structure. The SUp(3) breaking is taken into account by the strange quark mass 
which appears in the mass term, in the kinetic part and in the hyperfine part. The lat- 
ter also breaks SUp(3) through the masses of the pseudoscalar mesons exchanged among 
quarks. 
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2. The Hamiltonian 

The Hamiltonian has the form [ ITS] 

h = 5> + E Jr - f^+E^K) +E^;), (i) 

with the linear confining interaction 

VMj) = -~ K ■ \) C r y , (2) 
and the flavour-spin interaction 



^(%) = < £ V*(r«) Af Af 
If=i 

+ E V K(rij) Af Af + 7,(r y ) A?Aj + tyfa) A°A°) <?, • a, 
The analytic form of V 7 (r) (7 = ir,K,r] or 77') is 

with the parameters: 

flp = 9^ = sSca = o.67 §V = x 206, 

4-7T 47T 47T ' 47T 

r = 0.43 fm, a = 2.91 fin -1 ,*? = 0.474 /m" 2 , m M = 340 MeV, m s = 440 MeV, (5) 
fji n = 139 MeV, n v = 547 MeV. /v = 958 MeV, fi K = 495 MeV. 

which lead to a good description of low-energy non-strange and strange baryon spectra. 
Fixing the nucleon mass at mjv = 939 MeV, this parametrization gives for example tua 
= 1232 MeV and A/" (1440) = 1493 MeV. The lowest negative parity states appear at 
iV(1535) — iV(1520) = 1539 MeV, i. e. above the Roper resonance, in agreement with the 
experiment. 

3. The wave function Ansatz 

We start with the q subsystem and treat the quarks as identical particles in all cases. 
Then following Ref. [[12] the orbital (O) part of the lowest totally antisymmetric state 
must carry the symmetry [31] G . In the flavour-spin (FS) coupling scheme this state has 
the form 

|1> = |[31]o[211] c [llll] oc ; [22] F [22] 5 [4] F5 > (6) 

which means that the wave function is totally symmetric in the flavour-spin space and 
totally antisymmetric in the orbital-colour (OC) space and that the q 4 subsystem carries 
non-zero angular momentum and has zero spin. Then the q A q state is obtained by coupling 
the antiquark to the state |1) of Eq.© which leads to either 10^ or to 8^ and to a total 
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spin 1/2. To derive the orbital part we denote the quarks by 1, 2, 3 and 4 and the 
antiquark by 5 and introduce the internal Jacobi coordinates 
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(7) 



x = r i - f 2 , y = (fi + f 2 - 2f 3 ) / 
<? = (rl + f 2 + f 3 - 3f 4 ) /\/6 , t = (n + f 2 + r 3 + r 4 - 4f 5 ) / ^ 

The key issue is to construct a wave function with correct permutation symmetry in 
terms of the above Jacobi coordinates. Assuming an s 3 p structure for [31] G , the three 
independent [31] G states denoted by ipi are [ IT2] 
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e shell model wave functions and we took the quantum number m = 
everywhere, for convenience. Thus each function carries an angular momentum I = 1 in 
one of the relative coordinate, which leads to a total parity P = + 1 and a total angular 
momentum J = 1/2 or 3/2. The degeneracy of these two states can be lifted by the 
introduction of a spin-orbit coupling. 

The functions (Pl- lfTUj) are used to construct a totally antisymmetric orbital-colour 
state for the q 4 subsystem, in agreement with (jUJ). The coefficients of the resulting linear 
combination are fixed by group theory, namely by the Clebsch-Gordan coefficients of the 
permutation group S4. In this case, the absolute value of all three coefficients is equal to 
l/v^, which means that each of the states (fK|)- (fTT)|) contribute with equal probability. 

The pentaquark orbital wave functions are obtained by multiplying each by (t |000^> 
which describes the motion of the q A subsystem relative to q. The wave function associated 
to each relative coordinate is chosen to be a Gaussian. This gives 



^ = 
^2 = 
^3 = 
where 
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exp 
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(ii) 

(12) 
(13) 

(14) 



The two variational parameters are a, the same for all internal coordinates of the q 4 
subsystem, x, y or z, and (3, for t, the relative coordinate of q 4 to q. 

The algebraic structure of the state (JBJ) is identical to that of Ref . [ EH] ■ The small 
overlap of the resulting q 4 q state with the kinematically allowed final states could partly 
explain the narrowness of + . 
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Table 1 

The hyperfine interaction V x , Eq. Q, integrated in the flavour-spin space, for four 
quark subsystems. The upper index indicates the flavour of every interacting qq pair. 

g 4 I /, h I v* 

uudd 0,0 30 V v - 2 V™ - 4 VJ V 

uuds 1/2, 1/2 15 Kr - - 2 V^* + 12 + 2 - 2 V^ s 

ofofss 1, -1 V v + \V™ + |V^ U + + \V°, S + 2W K + fv v as - fV^ s 



4. Matrix elements 

The expectation values of the hyperfine interaction V x , Eq. (jSJ), in the flavour-spin 
space, are presented in Table Q for the three q 4 subsystems necessary to construct the 
antidecuplet. They are expressed in terms of the two-body radial form (j3J) now denoted 
as V^ aqb where q a qb specifies the flavour content of the interacting pair. The SU(3)^ 
is explicitly broken by the quark masses and by the meson masses. By taking VL MU = 
V^ s = V^ s and V^ u = V^, s = 0, one recovers the simpler model described in Ref. [ H§] 
where one does not distinguish between the uu, us or ss pairs in the ?7-meson exchange. 
Moreover, in Ref. [UU] one takes as parameters the already integrated two-body matrix 
elements of some radial part of the hyperfine interaction, as in Ref. Here we specify 

a radial form, which allows the explicit introduction of radial excitations at the quark 
level, whenever necessary. Then, from Table ^ one can easily reproduce Table 3 of [ E|] 
containing the coefficients x±,X2 and x^, i. e. the multiplicities, or the fraction of the two 
body matrix elements associated to 7T, K and r] exchange respectively, which appear in 
the expression for the mass. The first and last row of Xj, corresponding to G and H 
are straightforward, inasmuch as their contents are uudds and ddssu respectively. To get 
the Xi associated with N 5 and S 5 , which we call here Njq and S^- respectively, one must 
construct the linear combinations 

1 2 

^xC^io) = ^V x {uudd) + -V x {uuds), 

1 2 

in agreement with the flavour wave functions given in the Appendix and the relation 
V x {uuss) = V x (ddss). Moreover, in Ref. [H22], for each exchanged meson, one assumed 
that the radial two-body matrix elements are equal irrespective of the angular momentum 
of the state, I = or i = 1, which we won't do. 

5. Results and discussion 

In Table El we present the variational energy E of the model Hamiltonian (0) resulting 
from the trial wave function described by Eqs. (fTTj) -([Tl |) for various q 4 q systems related to 
the antidecuplet or the octet. One can see that, except for the confinement contribution 
(V c ), all the other terms break SU(3)jr, as expected: the mass term J2n=i m * increases, 
the kinetic energy (T) decreases and the short range attraction (V^.) decreases with the 
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Table 2 

Expectation values (MeV) and total energy E = ELi rra + (T) + (V c ) + (V x ) obtained 
from the Hamiltonian (JTJ) for various q A q systems. The mass M is obtained from E by 
subtraction of 510 MeV in order to fit the mass of B + . The values of the variational 
parameters a and j3 are indicated in the last two columns. 



4— 

q q 


ELi m i 


(T) 


(Vc) 


(v x ) 


E 


M 


at(fm) 


(3{fm) 


uuddd 


1700 


1864 


442 


-2044 


1962 


1452 


0.42 


0.92 


uudds 


1800 


1848 


461 


-2059 


2050 


1540 


0.42 


1.01 


uudsd 


1800 


1535 


461 


-1563 


2233 


1732 


0.45 


0.92 


uudss 


1900 


1634 


440 


-1663 


2310 


1800 


0.44 


0.87 


ddssu 


1900 


1418 


464 


-1310 


2472 


1962 


0.46 


0.92 


uusss 


2000 


1410 


452 


-1310 


2552 


2042 


0.46 


0.87 



Table 3 

The antidecuplet mass spectrum in MeV. 



Pentaquark 


Y, I, h 


Present result 


Carlson et al. [ITS'] 


Exp + GMO formula 


e+ 


2, 0, 


1540 


1540 


1540 




1, 1/2, 1/2 


1684 


1665 


1647 


^10 


0, 1, 1 


1829 


1786 


1755 




-1, 3/2, -3/2 


1962 


1906 


1862 



quark masses. For reasons explained in the introduction, we subtract 510 MeV from the 
total energy E in order to reproduce the experimental G + mass. 

For completeness, in the last two columns of Tabled we also indicate the values of the 
variational parameters a and (3 appearing in the radial wave function (fT4^1 which minimize 
the energy of the systems displayed in the first column. The parameter a takes values 
around «o = 0.44 fm. In the same quark model this is precisely the value which minimizes 
the ground state nucleon mass [E] when the trial wave function is <ft oc exp[— (x 2 +y 2 )/4al] 
where x and y are the first two of the Jacobi coordinates (J7J) defined above. The quantity 
«o gives a measure of the quark core size of the nucleon because it is its root-mean-square 
radius. The parameter j3 is related to the relative coordinate t between the center of mass 
of the q 4 subsystem and the antiquark. It takes values about twice larger than a, which 
is an indication that the four quarks cluster together, whereas q remains slightly separate 
in contrast to certain Ansaetze recently promulgated in the literature. 

Table El reproduces the calculated antidecuplet mass spectrum obtained from the mass 
M of Table El The masses of G + and H can be read off Table El directly. The other 
masses are obtained from the linear combinations 

1 2 

M(Njq) = -M(uuddd) + -M(uudss), 

o o 

2 1 

M(E W ) = -M(uudsd) + -M(uusss) . (16) 

3 3 
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In comparison with Carlson et al. where the mass of B + is also adjusted to the 

experimental value, we obtain somewhat higher masses for Njq, £jq and H , the latter 
being about 100 MeV above the experimentally found mass of 1862 MeV [[3j. This is 
in contrast to the strongly correlated diquark model of Jaffe and Wilczek [I2D]> where 
5 lies about 100 MeV below the experimental value. Note that the mass of + is also 
fixed in that model. In the lowest order of SU(3)^ breaking, one can parametrize the 
result by the Gell-Mann-Okubo (GMO) mass formula, M = Mjq + cY. In the present 
case one obtains M ~ 1829 — 145 Y. The fit to the measured masses of G + and H 
gives M ~ 1755 — 107 Y. Accordingly, the masses assigned to Nj^ and Ejq are 1647 
MeV and 1755 MeV. They are indicated in the last column of Table |3] Starting from 
this fit, Diakonov and Petrov [|^] analyzed the masses of the non-exotic members of the 
antidecuplet as a consequence of the octet-antidecuplet mixing due to SU(3)i? breaking. 2 
A new nucleon state at 1650-1690 MeV and a new £ at 1760-1810 MeV have been proposed 
as mainly antidecuplet baryons with Y = 1 and Y = respectively. Shortly after, Pakvasa 
and Suzuki [ 22] also considered the octet-antidecuplet mixing in a phenomenological way 
starting from the Gell-Mann-Okubo mass formulae. There, the resonance AT* (1710) was 
taken as the Y = 1 partner of G + , as in the original work of Ref. That analysis 

showed that the range of values for the mixing angle required by the mass spectrum of 
the Y = 1 baryons is not consistent with the range needed to fit strong decays. 3 

However the recent modified PWA analysis [|21j reconsiders the antidecuplet nature of 
AT* (1710) used to determine the mass of 6 + in Ref. [Ej. As a result, instead of AT* (1710), 
it proposes two narrow resonances 1680 MeV and/or 1730 MeV, as appropriate Y = 1 
partners of G + . This interpretation of the data clearly requires octet-antidecuplet mixing. 

In the present model, which contains SU(3)^ breaking, the mixing appears naturally 
and it can be derived dynamically starting from the Hamiltonian (0). Recall that Table 3, 
column 3 gives the pure antidecuplet masses. The pure octet masses are easily calculable 
using Tableland the octet wave functions (see Appendix). We obtain 

2 1 

M(N 8 ) = -M(uuddd) + -M(uudss) = 1568 MeV, 

3 3 

M(E 8 ) = -M(uudsd) + -M(uusss) = 1936 MeV. (17) 
3 3 

The octet-antidecuplet off-diagonal matrix element, denoted by V, has only two non- 
vanishing contributions, one coming from the mass (first) term of (0) and associated with 
the overlap of <&(Njq) and $(A r 8 ), or of $(£jq) and $(S 8 ), and the other coming from the 
hyperfine interaction. Using the Appendix one can obtain the analytic form of V as 

^{m s -m„) + | [Siuudss) - Siuuddd)] = 166 MeV for N 

(18) 

^(m s - m u ) + ^ [S(uusss) - S{uudsd)} = 155 MeV for £ 

2 A similar analysis, but restricted to the ideal mixing postulated by Jaffe and Wilczek [ , has been 
made in Ref. [123]. 

3 A more extended representation mixing including the 8, 10, 10, 27, 35 and 35 were considered in Ref. [ 
135) in the context of the chiral soliton model. The masses of Njq and Sjq were predicted to be the same 
as those in the last column of Table |3| The estimated range for the pure exotic pentaquarks turn out to 
be 1430 MeV < M(8+) < 1660 MeV and 1790 MeV < M(E~) < 1970 MeV. 
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where S = (T) + (V x ). The numerical values on the right hand side of Eq. (JTHJ) result 
from the quark masses given in Eqs. (j5J) and from the values of (T) and {V x ) exhibited in 
Table |21 One can see that the mass- induced breaking term is identical for N and E, as 
expected from simple SU(3) considerations. Its numerical value, 94.28 MeV, represents 
more than 1/2 of the total off-diagonal matrix element. 

The masses of the physical states, the "mainly octet" N* and the "mainly antidecuplet" 
N 5 , result from the diagonalization of a 2 x 2 matrix in each case. Accordingly the nucleon 
solutions are 

N* = N 8 cos 6 N - N w sin 6 N , 

N 5 = N 8 sm9 N + N w cos6 N , (19) 

with the mixing angle defined by 
2V 

tan 20 N = ——7-7 — r 7T7T7T • (20) 

The masses obtained from this mixing are 1451 MeV and 1801 MeV respectively and the 
mixing angle is 6^ = 35.34°, which means that the "mainly antidecuplet" state N 5 is 67 % 
Njq and 33 % N&, and the "mainly octet" N* the other way round. The latter is located 
in the Roper resonance mass region 1430 - 1470 MeV. However this is a q 4 q state, i. e. it 
is different from the g 3 radially excited state obtained in Ref. [E3 at 1493 MeV with the 
parameters (J3J) and assigned to the Roper resonance. A mixing of the q 3 and the q A q states 
could possibly be a better description of reality. There is some experimental evidence that 
two resonances, instead of one, separated by about 100 MeV, and located around 1440 
MeV, could consistently describe the tt — N and a — p scattering in this region [ 26J, 
however. Thus the issue of the existence of more than one resonance with J p = 1/2+ 111 
the 1430- 1500 MeV mass range remains unsettled. The "mainly antidecuplet" solution 
at 1801 MeV is far from the higher option of Ref. [121], at 1730 MeV, interpreted as the 
Y = 1 narrow resonance partner of + . 

In a similar way we obtain two £ resonances, the "mainly octet" one being at 1719 MeV 
and the "mainly antidecuplet" one at 2046 MeV. The octet-antidecuplet mixing angle is 
Oy, = —35.48°. The lower state is somewhat above the experimental mass range 1630 - 
1690 MeV of the the £(1660) resonance. As the higher mass region of £ is less known 
experimentally, it would be difficult to make an assignement for the higher state. 

The mixing angle 9n and 8-^ are nearly equal in absolute value, but they have opposite 
signs. The reason is that M(A%) > M(N 8 ) while M(E W ) < M(£ 8 ). Interestingly, each 
is close to the value of the ideal mixing angle 9 N = 35.26° and # s = —35.26°. Only the 
relative strengths of decays and selection rules can discriminate between mixing schemes 
as well as between models [ 231 12Z] ■ This is a task for a future work. 



6. Conclusions 

In conclusion we have used a variational method, which provides upper bounds on the 
masses of all isomultiplets of the pentaquark antidecuplet. We calculated dynamically the 
masses of the pure exotic pentaquarks + and 5 and the masses of the other members 
of the antidecuplet. The model on which these calculations are based reproduces well 
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the baryon spectrum, when baryons are described as q 3 systems. It assumes a flavour- 
spin structure for the hyperfine quark-quark interaction and its radial shape contains 
parameters which have been fitted not only to the ground state baryons, but also to a large 
number of excited states I n particular this interaction places the Roper resonance, 

modeled as a g 3 system, below the lowest negative parity baryons, in agreement with the 
experiment. However the description of strong decays in this model is not satisfactory 
(see e. g. Ref. [ 28J). Besides the qq interaction a qq interaction is necessary to describe 
pentaquarks. Here we did not introduce it explicitly but relied on the conclusion of Ref. [ 
ITB] that an attractive spin-spin interaction that operates only in the qq channel can lower 
the q 4 q energy to accommodate the + . In this way we can explain the mass shift of -510 
MeV necessary to reproduce the mass of + . It follows that this flavour- independent 
interaction equally lowers all the other members of the antidecuplet and of the octet. 

But in the new light shed by the pentaquark studies, the usual practice of hadron 
spectroscopy is expected to change. There are hints that the wave functions of some ex- 
cited states might contain q A q components. These components, if obtained quantitatively, 
would perhaps better explain the widths and mass shifts in the baryon resonances [ 29 . 
In particular the mass of the Roper resonance may be further shifted up or down. In that 
case the model parametrization should be revised and more precise four- and five-body 
calculations should be performed. On the other hand a full experimental confirmation of 
the G + and of the H resonances and more appropriate partial wave analysis of existing 
data would be of great help in understanding the structure of pentaquarks and of ordinary 
baryons. 



Appendix 

Here we give the form of one of the two independent flavour wave functions for each 
isomultiplet belonging to 10^. It is the function where both pairs of quarks, 12 and 34, are 
in an antisymmetric state (p[n](q a qb) = (q a qb — Iblal/V 7 ^- By analogy with the g 3 system, 
we shall use the label p for all states which are antisymmetric under the permutation (12). 
For G this wave function is straightforward 

$"(0) =<j ){11] {ud) ( t> [lx] {ud)s . (21) 

The Njq flavour wave function is obtained from that of by applying the [/-spin ladder 
operator LL of SU(3). Its normalized form becomes 

® P ( N To) = ^{[0[n]( MS )0[n]( MC O + 0[ii](«d)0[ii](ws)]s + <t>[n]{ud)(f)[n]{ud)d} . (22) 

Applying LL again one obtains the wave function of £y which is 
1 

$ p (S T q) = -j={(p [u] (us)(j)[n](us)s + [0[ii](ws)0[n] (wcf) + 0[u](ud)0[n](ws)]d} . (23) 

The wave function of H is as simple as that of but with another quark content of 
course 

$'(S— ) = [ i 1] (da)0 [ ii](ds)S ■ (24) 
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In these functions the normal order of particles 1234 is understood. In each case one can 
get the other linear independent function in the flavour space, $ A , with the quark pairs 12 
and 34 in a symmetric state, [2 ] (<?«<%) = (q a qb + q b qa)/^ ( Qa ^ Qb) or ^{qaQa) = QaQa, 
by applying the permutation (23) to the above corresponding function (see e. g. [ I5U|)- 
For example we have 

$ A (6) = J^[(f)[2]{uu)(j)[2]{dd) + <j)[2]{dd)<f)[2](uu) - (j>[2]{ud)<f)[2]{ud)\s . (25) 

Both the $ p and $ A functions are necessary in the calculation of the matrix elements of 
the hyperfme interaction. 

In the same notation, the N s and S 8 the flavour octet wave functions, antisymmetric 
under the permutation (12) are 

1 [2 
$ P (A 8 ) = -^=[0[n](«s)0[ii](wd) + <f)[u](ud)(f)[ii](us))s - y -<p[u](ud)(f)[u](ud)d . (26) 

[2 _ 1 

$ p ( s s) = y -(l)[u](us)(j) [n] (us)s - -j= [4> [n] (us)(f> [n] (ud) + <f>[u)(ud)^[n][us)]d . (27) 
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